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ABSTRACT 

A recent analysis of UV data from the Interface Region Imaging Spectrograph IRIS reports plasma 
“bombs” with temperatures near 8 x 10"^ K within the solar photosphere. This is a curious result, 
firstly because most bomb plasma pressures p (the largest reported case exceeds 10^ dyn cm“^) fall 
well below photospheric pressures (> 7 x 10^), and secondly, UV radiation cannot easily escape from 
the photosphere. In the present paper the IRIS data is independently analyzed. I find that the 
bombs arise from plasma originally at pressures between < 80 and 800 dyne cm“^ before explosion, 
i.e. between < 850 and 550 km above T500 = 1- This places the phenomenon’s origin in the low-mid 
chromosphere or above. I suggest that bomb spectra are more compatible with Alfvenic turbulence 
than with bi-directional reconnection jets. 

Subject headings: Sun: atmosphere 


1. INTRODUCTION 

iPeter et ^ (I20l4 henceforth P2014) recently reported 
on “Hot Explosions in the Cool Atmosphere of the Sun”. 
Throughout the paper the 8 x 10“^ K plasma is identi¬ 
fied as belonging to the solar photosphere. Two major 
pieces of evidence are presented by P20I4 to support the 
idea that the bombs originate from photospheric plasma. 
P2014 make quantitative estimates of the electron den¬ 
sity, leading them to conclude 

“A density analysis based on the IRIS ob¬ 
servations of O IV and Si IV shows that the 
bombs form in the photosphere” 

The reasoning is that such high density plasma pre-exists 
only in the deeper layers. The second argument con¬ 
cerns the presence of narrow absorption lines of Fe II 
and Ni II seen against the bright emission. Cool plasma 
exists above the emitting plasma with sufficient opacity 
to produce absorption lines. P20I4 did not analyze these 
lines quantitatively, this is done below. 

In this paper I demonstrate that the Ri 8 x 10"^ K radi¬ 
ating plasma must exist at least 550 km above the contin¬ 
uum photosphere, above the canonical temperature min¬ 
imum region. I take the essential observations of P2014 
and re-analyze them within the following simple picture. 
I adopt, with justification, a ID model of the atmosphere 
that captures the essential stratification of 3D dynami¬ 
cal atmospheres to study the transfer of UV continuum 
and line radiation, in both the emission (C II, Si IV, 
Mg II) and absorption features (Fe II and Ni II) reported 
by P20I4. I assume that electrons are fully thermalized 
owing to the high densities encountered. I re-analyze 
the physics determining the ratios of permitted Si IV 
and spin-forbidden O IV] lines (“]” denotes, by conven¬ 
tion, a spin-forbidden transition), and the implications 
of the absorption features. I assume bomb spectra origi¬ 
nate from pre-existing plasma that is heated to produce 
the 8 X 10“^ K plasma, keeping the density approximately 
constant. The latter is assumed to be due to strong mag¬ 
netic forces, to be consistent with the 5 minute lifetimes 
of the bombs that greatly exceed dynamical time scales 


(r^ 1000 km / 100 km s ^). Finally I review the analysis 
of P2014. 

2. ANALYSIS 

Below I use phot oionization cr oss sections from the 
OPACITY project I SeatonI Il987tl . other atomic data 
compiled bv lJudgj ( 2007I1. IRIS data from P 2014 and 
calibration data from iDe Pontieu et al.l (I2014D . Photo¬ 
spheric abundances from lAllenI (|1973D are adopted. 

2.1. Pre-existing stratification and UV opacities 

Under undisturbed conditions, the solar photosphere 
and chromosphere are strongly stratified. I therefore dis¬ 
cuss bombs in the context of an average atmospheri c 
model, model C (“VAL3C”) of IVernazza et al.l (1198111 . 
This is a static ID model. Nevertheless it captures typ¬ 
ical density and pressure stratification of dynamical 3D 
models. Much has recently been made of the essential 
dyna mical nature of the chromosphere (jLeenaarts et ahl 
I2OI2I e.g.), but such dynamical departures from hydro¬ 
static equilibrium are large only in the tenuous regions 
of the upper chromosphere. Figure [1] shows the run of 
pressure with height z of the VAL3C model together with 
a probability distribution function derived from the 3D 
Bifrost model “en024048_hion” developed for the IRIS 
projeclQ. It is clear that VAL3C can be used to study 
the vertical transport of radiation, under undisturbed 
conditions, low in the solar chromosphere. 

The atmospheres of G-type stars contain enormous 
opacities from abundant neutral atoms: Rayleigh scat¬ 
tering from hydrog en and bou nd-free transitions from 
Si, C, Fe, A1 (e.g. lAllenI 1197^ . UV radiation cannot 
escape from or penetrate the undisturbed photosphere. 
Below 152 nm, the opacity in the lower chromosphere is 
dominated by bound-free opacity from the ground levels 
(3p^ ^Tb.1^2) of neutral silicon. The Si I photoionization 
cross section is crsi = 6 x 10“^^ cm“^. Using initially 
''^Si i/nsi ~ 1 (silicon is mostly neutral), and an abun- 

^ http://sdc.uio.no/search/simulations 
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Figure 1. The stratification of the VAL3C model is shown to¬ 
gether with the average and PDF of the stratifications from one 
random time step (7^491) taken from the Bifrost simulation. The 
lower-middle chromosphere is extremely well represented by the 
VAL3C stratification up to the point near 2000 km where dynam¬ 
ical effects start to dominate. Above this height departures from 
hydrostatic equilibrium are especially large and have drawn much 
attention in the literature. 


Table 1 

Optical depths and escape fractions for model VAL3C 
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6Xp( 4nnm ) 
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dyne cm“^ 
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0.74 
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2 X 10“® 

100 

1 
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725 
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dance nsi/’^H = 3.5 x 10 we find 


k = riH-crsi 
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»H 

10 ^^ cm“3 


cm 


( 1 ) 


With these opacities, optical depths r and the escape 
probabilities exp(—r) for line radiation lying in the Si I 
continuum are listed in Table 1. Below 600 km, where 
gas pressures p exceed 600 dyne cm“^, the optical depth 
(proportional to column mass m where p = mg in hy¬ 
drostatic conditions and g is the solar acceleration due 
to gravity), fewer than 1 in 100 photons escape from 
this high density, opaque plasma, dropping as a double¬ 
exponential function of height. 

To construct Table 1 I assume unrealistically that 
bomb radiation does not affect the plasma surrounding 
it. Photoionization of neutrals reduces UV opacity at 
140nm, but this effect is offset partly by increased re¬ 
combination as more electrons are made free. 

In principle a non-LTE transfer calculation is needed to 
investigate the statistical equilibrium of the bomb and its 
effects on the surrounding atmosphere. However, bomb 
plasma with densities > 10^^ cm“^ in a ID atmosphere 
is so geometrically thin (< 0.1 km) and hot that cal- 
culations with the RH non-LTE code (|Uitenbroekll200(1l) 
proved impossible to make converge. 

To proceed, two simpler calculations were made. First, 
RH was used with downgoing radiative intensities of 
transition region lines that equal the brightest observed 


bomb intensities (as shown below), as an upper boundary 
condition. The non-LTE rate equations were solved for 
H, C, O, Si, Fe to study the penetration of such bright ra¬ 
diation downwards and changes in ionization balance and 
the electron number densities, keeping the temperature 
fixed. Less than 10% of the 139.3 nm radiation pene¬ 
trates to or escapes from 2 ; = 500 km. To understand 
this result, I then made the following rough calculations, 
examining the consequence of placing the bright bomb 
emission immediately below the temperature minimum 
region (z = 450 km, p(z) = 1600 dyne cm“^). The undis¬ 
turbed continuum optical depth at 140nm is 100. The 
wavelength-integrated intensity of radiation in bomb B-1 
in the Si IV lines is Ib-i ~ 2.6 x 10® erg cm“^ s“^ sr“^. 
The photoionization rate of Si I by the Si IV bomb emis¬ 
sion is 

Pi, = 27T^a « 70 s-i (2) 

riv 

At z = 450 km the VAL3C electron density is ~ 4 x 10^^ 
cm“®, and the recombination rate per electron a ~ 2.5 x 
10“^^ cm®s“^. Ionization equilibrium gives 

^■Si i/«Si II = aUe/PiK « 1/70, (3) 

i.e. the 140nm continuum optical depth is reduced to 
100/70=1.4. The observed C II lines are almost as strong 
as the Si IV lines, and H La will certainly be strong. But 
all additional UV radiation will increase the electron den¬ 
sity as elements become ionized: fully ionized silicon and 
carbon contribute an additional 1.4 x 10^^ and 10^^ elec¬ 
trons per cm® at z = 450 km. With these additional 
electrons, equation ([3|) yields nsi i/nsi ii ~ 1/16 with 
a continuum optical depth of about 6, stopping all but 
2.5% of the radiation from emerging. These calculations 
explain the RH non-LTE calculations, so that to con¬ 
clude, UV bomb emission below 450 km is not visible to 
IRIS even in the most extreme case B-1. 

2.2. Density analysis 

When the radiation in two emission lines depends dif¬ 
ferently on electron density Ue and temperature Te, it 
is obvious that one line ratio can only broadly con¬ 
strain possible mean densities in the emitting plasma. 
The general problem of finding electron densities from 
thermal but low d ensity plasmas has been discussed by 
lJudge et al.l (I1997D . Refer to this paper for the notation 
used here. Under standard assumptions each line £ emits 
a total intensity of 

n oo 

fj,{ne,Te)Hi{ne,Te) dTedug erg cm“^ s“® sr“\ 

( 4 ) 

where /i is a differential emission measure and Hi de¬ 
pends on atomic processes. To derive a mean density n 
one must specify a form of the source term p, for ex¬ 
ample in an isothermal plasma at temperature Tq, p = 
aS{ne — n)S{Te — Tg). Lines of Si IV and O IV] have dif¬ 
ferent peak ionization equilibrium temperature tempera¬ 
tures, Tsi and Tq. P2014 look for solutions for n at con¬ 
stant pressure nTe =constant, so that p(n,Te) = p{Te) 
and where 

pin, Te) = pin) = aSin - Tsi) + bSin - Tq) (5) 
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Figure 2. The figure compares the ionization balance of Si and 
O ions in response to changes in the electron density. In the “coro¬ 
nal approximation” (low densities) the curves do not vary with 
electron density. However at densities higher than 10® cm“® the 
effects of metastable level populations and reduced dielectronic re¬ 
combination rates significantly change, systematically, the relative 
ion populations. Notably, the Na-like Si IV ion’s population grows 
rapidly with increasing density. The opposite is true for the B-like 
O IV ion. 

Note that the line ratio yields the mean electron density n 
only when the ratio of the emission measures a at Tsi and 
b at To is known. In this regard the calculations of P2014 
are not specified sufficiently to reproduce them. I find 
them to be consistent with a/6 ~ 4 using calculations 
with similar ionization equilibria. It appears then that 
they assume the geometrical thickness of the emitting 
regions is the same in the Si and O lines, so that the Si 
line has roughly 4x the emission measure (i.e. a/6 = 4 
in the present notation). 

The specific case for Si IV/0 IV] intensity ratios has 
additional problems: 

1. Si is a low-first ionization potential (low-FIP) ele¬ 
ment and O is not. Differential abundances Aqi/Ao 
between low- and high FIP e lements are well doc¬ 
umented fe.g. lFeldmaijlI99^ . Si can be enhanced 
by factors of up to 3 compared to O. 

2. The Si3+ and 0^+ ions have different energy 
thresholds for ionization from the neutral species 
(58.7 and 103.7 eV). Consequently ionization equi¬ 
librium temperatures differ by a factor of two 
(P2014 use logipTe = 4.87 and 5.16). 

3. Ion populations depend on micro-fields in finite 
density plasmas (above ~ 10® cm“®, iSummeri 
[l97l . Above 10® cm ® collisions reduce rates 
of dielectronic recombination (DR) by perturbing 
the doubly excited sta tes before the A uger stabiliz¬ 
ing transition occurs (|Summerslll974f) . In Na-like 
Si®"*", DR to Si^+ is significantly suppressed (see 
the computations of the homologous Li-like ions by 
IDovle et al.ll200^ . Together with increased ioniza¬ 
tion from metastable levels of B- and Be-like ions, 
at plasma densities above 10®^ cm ~® I estimate us¬ 
ing the approach of lJuded (|2007ll that the emitted 
power of Na-like Si IV lines above those of O IV] 
lines by a factor of three to four (Figure [5]). 

4. The Si IV line core intensity ratios 139.3/140.2 are 
close to but (at least for bomb B-1 where the ratio 
is 1.6-1.7)) less than 2:1. Scattering in the out¬ 
ward (observed) direction increa ses the 139.3 com- 
ponent’s integrated intensity (|Judge fc Pietaril^ 


1200411 since the lines are not optically thin, unlike 
the O IV] lines. In essence, downward directed ra¬ 
diation in Si IV lines can be scattered outwards, 
but in O IV] lines downward radiation is absorbed 
by the background continuum. 

Another problem is that intensities of UV lines of Na- 
like and Li-like ions are known to be too weak by factors 
of several compared with lines of other sequences. I n the 
accurate measurements used by l|Judge et al.lll995l) . the 
Na-like Si IV lines are 2.5 times brighter than calcula¬ 
tions predict. This discrepancy must lie outside of the 
high density effects because disk-averaged spectra con¬ 
tain many spin forbidden lines whose intensities are in¬ 
consistent with rig > 10®° cm“®. Perhaps this originates 
in non-equilibrium ionization (iHansteenl 1199311 or non- 
thermal distributions of electrons. Here I simply remind 
the reader of this additional problem. 

Taken together, all these issues indicate that compu¬ 
tations of relative intensities of Si IV and O IV] lines 
are uncertain to at least a factor of three. Most impor¬ 
tantly, at least four independent processes can enhance 
the intensities of the Si IV line, each by factors of two to 
three. When these systematic effects are compounded, 
low-density computations may underestimate the ratio 
of Si IV to O IV] line intensities by an order of magni¬ 
tude, which then would be mistakenly interpreted as an 
increased electron density using the above methods. 

Figure [3] (a) shows Figure S8 of P2014, generalized 
to include different assumed values of a/6. Differences 
with the shape of the curve of P2014 are primarily due 
to the different ionization equilibrium calculations used. 
I included app roxima t e corr ections for the DR suppres¬ 
sion following lJudgel (120071) . Clearly, the chosen form 
of ^i{Tg) strongly determines the outeome of the density 
analysis. There is no physical reason for choosing par¬ 
ticular values of a/6 in the bombs. When electron heat 
conduction from overlying coronal plasma dominates the 
energy balan ce, these quantities are physically related 
(|Jordanlll992D . But there is no evidence for coronal emis¬ 
sion in the bombs, and there exists a class of solutions 
for “cool loop” structures invoked to account for long- 
standin g observational problems in the so lar transition 
region (iFeldmanl 119831: iDowdv et al.lll986D . Such solu¬ 
tions can in principle have ratios a/6 that greatly exceed 
unity. 

Figure [3] (b) shows line ratios computed for homoge¬ 
neous plasmas where ^(ug, Tg) = 6(ng — no, Tg — Tq). All 
regions below the contour marked with the arrow are 
compatible with the observed lower limit to the ratio. 
When 6 <C a (very little plasma above Tg = 10® K) the 
line ratios contours become more horizontal than verti¬ 
cal. The ratios are then primarily sensitive to Tg than ng. 
In adopting equal emission measures for T = 4.87 and 
5.16, P2014 are forced to look for solutions beyond 10®® 
cm“® in their Figure S8 in order to reconcile observed and 
computed line ratios. Instead, Figure [3] demonstrates 
that the line ratio analysis is subject to an entirely dif¬ 
ferent interpretation where there is less plasma above 
logioT/ = 4.9. 

In conclusion, the mean density is not even defined 
without drastic simplifying assumptions, even then they 
are poorly constrained and subject to systematic over¬ 
estimation by factors of 3 if low density ionization equi- 
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Figure 3. The computed ratios of line emission is plotted as a 
function of electron density Tie for the lines Si IV 139.3 nm to 
O IV 140.1 nm, assuming constant pressure. The densities for the 
temperature of the Si IV line are shown. The ratios are computed 
with Equations [4| and [5] A ratio of 200, a typical observed value 
for the bombs, is shown as a horizontal line with an arrow, (b) 
Computed ratios are plotted as a function of electron temperature 
Te and density rig. Homogeneous plasmas were assumed. Lines of 
Si IV amd O IV form at different temperatures so that there is a 
gradient in the ratio dRjdT below logj^Q = 5.1. 


Table 2 
Absorption lines 


Ion 

A 

nm 

log gf 

^lower 

cm~^ 

P(to = 1) 
dyn cm“^ 

S I 

139.2589 

-3.35 

573 

360 


140.1515 

-1.18 

396 

2.3 

Fe II 

133.5409 

-0.93 

27,620 

800 


139.2816* 

-0.483 

20,806 

60 


139.3214 

-1.21 

20,340 

290 


140.1774 

-0.91 

22,810 

250 


140.3101* 

-1.922 

1,872 

25 


140.3255 

-1.25 

22,810 

550 

Ni II 

133.5203* 

-0.19 

1,507 

4 


139.333* 

-0.88 

0 
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requires the highest optical depths to form a weak ab¬ 
sorption line, the pressure at tq = 1 being almost 800 
dyne cm“^. This corresponds to a column mass m of 0.03 
g cm“^ or a height of about 560 km. The line probably 
requires an optical depth close to 0.1 or even less since it 
is very shallow (bomb B-1) or absent (other bombs), thus 
these absorption features must form in plasma originally 
above about 80 dyne cm“^, near a pre-existing height 
of 850 km. The absorption of each line is stronger for 
smaller values of p(ro = 1), The relative weakness of the 
observed S I lines implies that S I is mostly ionized above 
this height, a result confirmed in radiative transfer cal¬ 
culations with RH. Other than this, the line strengths 
appear in qualitative agreement with the absorption line 
depths seen in the data of P2014 (their Figure 3). 

3. DISCUSSION 

I have summarized and/or demonstrated the following: 

• UV photons from Si IV bombs cannot escape from 
below 500 km above the continuum photosphere. 

• Computed Si IV lines have been shown to under¬ 
estimate their brightness by factors of 2-3 in much 
earlier work, even in quiet Sun plasmas. 

• At very high densities claimed by P2014, Si IV be¬ 
comes yet stronger as dielectronic recombination to 
Si2+ is suppressed and ionization from metastable 
levels increases. Such effects were not taken into 
account by P2014. 

• In at least one case the Si IV lines show evidence 
of a small optical depth. Scattering of Si IV pho¬ 
tons by modestly opaque bomb plasma will en- 
han ce outward directed inte nsity by up to a factor 
of 2 lJud£fe fc Pietaril^l2004 


*Data from Kurucz CD-ROM 23. Other data from NIST’s spec¬ 
troscopic database. Partition functions of 50 and 25 were used for 
Fe II and Ni II corresponding to a plasma temperature of 6500K. 

librium approximations are used. By acknowledging pos¬ 
sible FIP effects and the anomalous Li- and Na- like ion 
intensities, the over-estimate can easily increase to an or¬ 
der of magnitude or more as the increased Si/0 intensity 
ratio is ascribed erroneously to pure density increases 
that decrease the O IV] emission rate. 

2.3. Absorption lines 

P2014 report several lines seen weakly in absorption 
against the bomb emission. They find that they are 
all slightly blue-shifted owing to flux emergence push¬ 
ing pre-existing material upwards, and conclude that 
“the absorption features indeed belong directly to the 
bombs”. Here I simply estimate the column mass of ma¬ 
terial needed to account for these weak absorption fea¬ 
tures. 

Table 2 lists pressures in the VAL3C atmosphere where 
the optical depths in the background atmosphere equal 
one in those observed absorption lines for which NIST 
atomic data are available. (Optical depth is proportional 
to column mass m and p = mg). Both Fe II and S I 
were assumed to be dominant ionization stages in the 
chromosphere, in making this table. The 133.54 nm line 


• Any density analysis of Si/0 line ratios hinges on 
an arbitrary assumption concerning the form of the 
emission measure distribution between log;^Q Tg = 
4.8 and 5.1. Unfortunately, almost any density can 
be made consistent with the data, given this extra 
degree of freedom. 

• IRIS should observe the permitted multiplet of 
O IV near 134 nm if it is to address this problem. 

• Absorption lines with optical depths O.I form in 
normal chromospheric plasma at pressures close to 
80 dyne cm“^ which originally lies near 850 km. 

Let us give P2014 the benefit of the doubt and ac¬ 
cept their general assumptions in evaluating densities. 
Even if the assumptions concerning the form of /r(ne, Te) 
are accepted, their calculations miss important physi¬ 
cal processes leading them to under-estimate the inten¬ 
sities of the Si IV lines, by a factor of 10 or so. Conse¬ 
quently they over-estimate the electron densities by an 
order of magnitude. Even in the most extreme case (B- 
I), log;^Q rig ^ 12.7. This new constraint implies most 
bomb plasma originated from regions where, before ex¬ 
plosion, p > 10“^^ g cm“^ or a height < 1180 km above 
the undisturbed continuum photosphere. This constraint 
may seem uncomfortably above the 850 km limit derived 
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from absorption lines in bomb B-1. But given the general 
difficulties with density diagnosis, and the fact that the 
850 km limit comes from the weakest of the absorption 
lines visible only in B-1, it seems reasonable to conclude 
that, all constraints considered, the bombs originate from 
plasma between 500 and > 850 km in the pre-bomb atmo¬ 
sphere. If their implicit assumption {a/h = 4) proves in¬ 
correct, then the electron and hydrogen densities will be 
change accordingly (Figure |3ljb)). No method is known 
to the present author by which this important assump¬ 
tion is testable given these data. New observations in¬ 
cluding the permitted O IV multiplet near 134 nm should 
help. 

Are the lower limit pressure and height of 500 km 
“photospheric?” Certainly it lies above the temperature 
minimum in VAL3C. However, as understood in physical 
terms, the photosphere is the region where the bulk of 
the radiative flux is radiated into space. The solar photo¬ 
sphere starts therefore where the radial optical depth in 
the 500 nm continuum equals unity, T 500 = 1. The upper 
levels of the photosphere are less well defined. Allen’s 
(1973, § 77) definition is that the photosphere extends 
to where rsoo = 0.005 which lies 320 km above T 500 = 1. 
Only a small fraction T 500 < 0.005 of the Sun’s radiative 
flux can originate above tsqo = 0.005, so Allen’s defini¬ 
tion allows > 99.5% of the solar flux to be radiated. 

Does it even matter that the bombs form in “pho¬ 
tospheric” versus “chromospheric” plasma? Or is this 
merely semanticism? Indeed it matters because of the 
stratification. At the high densities of the photosphere, 
the Alfven speed is slow for field strengths of a few hun¬ 
dred Gauss, characteristic of the regions without spots or 
pores where these bombs appear. Using Fig 1 of P2014, 
I use (optimistically) 5=400 G. The Alfven speed Va at 
500 km is 9 km s“^. Fast reconnection tends to occur 
at rates with outflows close to O.IVa. For a bomb of 
say 3Mm size, the reconnection time scale is of order an 
hour, an order of magnitude larger than a typical bomb 
lifetime. 

All things considered, it seems that the bombs arise 
from plasma with total particle densities, before explo¬ 
sion, that lie between 550 and > 850 km above the 
T'soo = 1- This places the phenomenon’s origin in the 
low-mid chromosphere or above. 

4. SPECULATIONS 

Prompted by the higher Alfven speeds associated with 
chromospheric plasma, I speculate on the nature of the 
bomb plasma. I also note that other observational prop¬ 
erties of the bombs identified by P2014 are not easy to 
reconcile with the association of bombs with reconnec¬ 
tion jets. More generally, reconnection leads to efficient 
acceleration of plasma through the change in topology 
and subsequent Lorentz force, but the direct heating ef¬ 
fects are usually small. Why then should we expect to 
interpret the radiation from strongly heated plasma such 
as seen in the IRIS emission lines directly as reconnection 


jets? 

The 5 minute bomb lifetimes greatly exceed the dy¬ 
namical crossing times if we are to accept the proposal 
that the 100 km/s speeds represent reconnection jets. 
Reconnection must occur over several minutes in this pic¬ 
ture, a perfectly reasonable picture. However, the bombs 
would be expected to exhibit large asymmetries between 
up- and down-flowing plasma because the surrounding 
atmosphere is highly stratified. The proposed supersonic 
jets must surely interact with surrounding magnetized 
plasma with observable consequences. At least for the 
overlying chromosphere, no obvious interaction is seen, 
instead benign absorption features are observed. What 
happens to the energy of the proposed jets? Perhaps this 
energy returns to photospheric layers, guided by closed 
fields, where it is deposited but is almost invisible be¬ 
ing energetically a small contribution. I suggest that a 
more natural explanation of the broad rather symmetric 
lines is Alfvenic turbulence trapped in a compact loop¬ 
like flux tube that exists somewhere in the low-middle 
chromosphere, driven perhaps by reconnection that oc¬ 
cupies a small volume of the system at any given time. 
Bombs are usually associated with magnetic footpoints of 
opposite polarities separated by 1000-2000 km, suggest¬ 
ing such a structure. This speculation appears worthy of 
further study. 
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